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INTRODUCTION 


HARMONIC ANALYZER is a frequency spectrum analysis package designed with the 
engineer in mind. It applies the concept of the Fast Fourier Transform (FFT) to 
an input data set to provide a frequency domain representation of the function 
approximated by that input data. | HARMONIC ANALYZER is based on DYNACOMP's 
popular FOURIER ANALYZER, but includes special data handling features. 


The user may save and recall both the input data file and previously calculated 
results. Data files may be loaded (from cassette, diskette, or disk), added to, 
deleted from, and generally edited. As with FOURIER ANALYZER, the input data may 
be plotted before the calculations are started. The frequency domain results may 
also be plotted as well as saved (to be later recalled for new plots). The 
analyses are easily repeated without re-entering the data. This is handy when 
you wish to test the sensitivity of the results to a particular data point. 


One of the unique features of HARMONIC ANALYZER is that the data need not be 
equally spaced in "time." A special spline interpolation is performed just prior 
to the FFT calculation to convert the the unequally spaced data into a new "data" 
file which is uniformly spaced as required by the FFT. 


RELATION TO OTHER DYNACOMP FREQUENCY TRANSFORM PROGRAMS 


DYNACOMP's FOURIER ANALYZER and TRANSFER FUNCTION ANALYZER were designed to be 
used on limited duration signals. The resulting transforms, therefore, represent 
a continuous spectrum. HARMONIC ANALYZER is meant to handle another class of 
signals; those which are repetitive and which, therefore, have discontinuous 
spectra. If FOURIER ANALYZER or TFA are used on single period waveforms 
belonging to repetitive wave-trains, the results will be wrong (not the fault of 
the program). Similarly, if HARMONIC ANALYZER is applied to limited duration 
signals, the results are very misleading. The truncated sinewave example given 
later demonstrates this. 


THEORY 


Repetitive signals may be decomposed into sinusoidal components as follows: 
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where 


Wy = ny (radians/second) 
Wy = 21/7, (radians/second) 
Te = period of signal (seconds) 


The Fourier Transform is defined as: 
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F(W) = Jeter dt 


oo 


Inserting the series representing f(t), we get: 
F(W) = an (W=W,) 
= 0 (otherwise) 
The Fourier Transform of repetitive signals is, therefore, discontinuous; com- 
posed of spikes. 
As a theoretical example, consider the following function: 


f(t) = 43 -1 st $1 


f(t) = 0 elsewhere 
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This is not a repetitive waveform. 


The associated Fourier Transform is a “sinc" function: 


F(w) = sin(w) 


Ww 
Note that this is a continuous function. 
If instead the function were: 
F(t) = % -1ét $1 
repeating at t=0, t=nT, 


> 
where ie = 2 


The spectrum would be discrete: 


e _ sin(W) 
F, = F(W,) = n 
Wh 


The larger Ty is, the closer together will be the Whe and in the limit as T,?@the 
spectrum will become continuous. 


A classical example of the difference between the Fourier integral for limited 
duration signals and the DFT/FFT is the spectrum derived for a limited duration 
sinewave. The DFT and FFT algorithms generally give a spectrum with a simple 
spike at the sinewave frequency. This is what one would expect for a continuous 
(infinite in extent) sinewave. The Fourier integral instead presents a spectrum 
in which the signal energy is concentrated about the sinewave frequency, but with 
some spread (actually an offset "sinc" function). 
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APPROXIMATION ERROR 


By using only a finite number of data points to represent a (piece-wise) 
continuous function, we are introducing approximation error. In effect, the 
function is being digitized*, or sampled. Thus, it is not possible to accurately 
and unambiguously capture the high frequency components of the original function. 
This limitation is described by the Nyquist sampling theorem which, in essence, 
states that the sampling rate (number of samples/unit time) should be at least 
twice the highest frequency component contained in the signal. For sampling 
rates below this, "aliasing" can occur. Spectrum calculations beyond this point 
should be considered suspect, though not necessarily invalid. Usually one avoids 
this problem by "filtering" the data before sampling. This may be done using the 
DATA SMOOTHER program also available through DYNACOMP, or DIGITAL FILTER (avail- 
able February 1981). If the data input to HARMONIC ANALYZER is equally spaced, 
the sampling rate will be twice the highest frequency calculated, though perhaps 
not at or greater than twice the rate of the highest component in the original 
signal. 


THE NUMBER OF DATA POINTS 


The FFT algorithm employed in HARMONIC ANALYZER requires that the number of data 
points processed be equal to an even power of 2 (e.g., 2, 4, 8, 16, 32, etc.). 
You will not be permitted to enter the transform step unless that is the case. 
The data file structure used by HARMONIC ANALYZER is generally the same as that 
employed by TFA (Transfer Function Analyzer). Thus, it is possible to load a 
data set which violates this condition, but data will have to be added or deleted 
to bring the number of points to an even power of two. 


NON-UNIFORM, DISORDERED DATA 


Before HARMONIC ANALYZER enters the transform stage, a new data set is con- 
structed from the one supplied. First, the "old" data set is ordered in terms of 
ascending "X" (or "t"). Next, the range of the N data points is determined and a 
new set of N equally spaced "X" (or "t") values chosen. Using a cubic spline 
interpolation, the "Y" (of "f") values corresponding to these equally spaced, 


* Actually, digitization is a further source of error. 
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ordered values are then determined. This artificial data set is subsequently 
sent to the FFT subroutine. If the original data set was ordered and equally 
spaced, the artificial data set is just an image of the original data collection. 


A cubic spline interpolation scheme was chosen over a cubic LaGrange (or Newton) 
in order to avoid the "ringing" associated with the latter. 


USING HARMONIC ANALYZER TO CALCULATE FREQUENCY RESPONSE 


HARMONIC ANALYZER may be directly applied to analyze the transfer function of a 
linear system. The experimental component is simple. A repetitive string of 
delta functions (e.g., pulse train of "spikes") is used as the input to the 
system to be analyzed. If the spike is transmitted unaltered, then F(w.) =a 
constant. The Fourier transform of the delta function is unity. Therefore, if 
the spike is transmitted with distortion, the transform (as calculated by 
HARMONIC ANALYZER) of this distorted waveform is proportional to the frequency 
response of the system being tested. An example is given in the attached 
listings. 


A word of caution is in order. Be careful that the input pulse is not so large 
that non-linear distortion occurs, or too wide. Ideally, the width should be 
less than the distance between sampling points. 


PLOTTING WIDTH 


The aspect ratio of the plotting routine is controlled by the parameter U. By 
increasing the value of U, fewer line feeds are inserted between plotted points. 
U is defined at the very beginning of the program and may be adjusted by the user 
to suit the terminal device employed. 


In conclusion, HARMONIC ANALYZER is a fairly general and very powerful frequency 
analysis tool. To get the best results from its application, remember that it is 
a numerical approximation whose accuracy and resolution improves with more 
precise and extensive data. 


HARMONIC ANALY: AMP 


The following set of examples was generated using the 
North Star version of HARMONIC ANALYZER. All other versions 


are very similar in thelr behavior. 
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THIS PROGRAM TRANSFORMS A GIVEN 
REPETITIVE DATA SET INTO ITS 
FREQUENCY DOMAIN REPRESENTATION 
USING THE FAST FOURIER TRANSFORM. 


THE USER INPUTS THE NUMBER 

OF DATA POINTS TO BE TREATED. 
THIS NUMBER MUST BE AN EVEN POWER 
OF TWO. FOR EXAMPLE? N=27K, 

THE PROGRAM WILL THEN PLOT 

THE DATA AND ITS SPECTRUM IN 
TERMS OF AMPLITUDE AND PHASE. 


YOU MAY ALSO CHOOSE TO RECALL 
PREVIOUS DATA OR RESULTS. 


DO YOU WISH TO PLOT 
PREVIOUS RESULTS (V/N)t 7M 


INPUT NUMBER OF TABLE VALUES! 716 


IS THE DATA TO BE LOADED 
FROM DISK (Y/N) 7N 


INPUT THE DATA. YOU MAY EDIT IT LATER. 


xOn 

YO 1) = 79 
xC2) = 12 

YO 2) = 79 
x3) = 72 

YO 3) = 19 
XC 4) = Th 

YO 4) = 70 
XC5) = 75 

YO S) = 79 
xe 6) = 76 

YO 6) = Th 
XC 7) = 17 

vO?) = 7h 
x @) = 7g 

yoo) = 71 
XC9) = 79 

Yoo) = th 
XC 10) = 740 

YC 100 = TL 
XC ADD = TBE 

yO 19) = 1h 
XC120 = The 

YO 12) = 
XC 13) = 745, 

YC 13) = 79 
XC 14) = TL6 

YC 14) = 79 
XC 45) = 747 

YO 13) = 70 
XC 16) = 712 


YC 16) = 79 


E.G., 2, &, 8, 16, 32 


The phase Is restricted to the range 0-360" 


A“Y" answer will lead to prompts concerning 
loading @ previously saved "results" file. 
Control will then go directly to the plotting 
routines for amplitude and phase. 


You will not be able to pass this point with 
an I legal number of data pairs. 


If the dete were to be loaded from disk (or 
cassette), the number of table values chosen 
above must be larger than the number of pairs 
In the file to be loaded. In this case we will 
build # 16 point data table and save It later, 


At this stage data Is entered unt!! 16 pairs ere 
recorded. If a mistake Is mede It can be corrected 
later. Note thet the date need not be entered In 
“time order, It will be sorted later. 


The example we will first examine Is a simple 
train of broad pulses. 


This Is an Intentional error to be corrected. 


Observe that this data polnt Is out of order. 


° 
I o---1 
* 1 -.375 
t t 
* 2 2 ~6375 
t t 
7 : 3 -.375 
' t 
* a) -.375 
! t 
* ts 6375 The ordered data set Is plotted 
$ t without modification. The values 
t * t 6 +625 printed to the right ere the result 
: 3 of removing the DC component. The 
t ® +7 +625 printed values are the ones sent to the 
: t Interpolation routine. 
* a) +625 
t t 
t * 29 +625 
t a 
t * 2 10 +625 Always try to use a plot width value 
: 1 greater than the number of data pairs. 
t ® ie +625 otherwise plot velues may have to 
t t double-up. If more than one point Is 
2 t 12 -.378 plotted on a given line, it Is not 
: t possible to also display the values 
: i413 -.375 to the right. 
* t 14 -.375 
sy : A maximum of two plot values will be 
* 2 15 -+375 shown on a given plot line. The third 
u : - a5 and higher will be ignored, 
I----1----I----I----I----I----1 
CONTINUE?Y, 
FOURIER TRANSFORM UNDER WAY.... 
FREQUENCY AMPLITUDE PHASE 
CYCLES/UNIT SCALED DEGREES 
° ° There Is somewhat of @ pause before the values are 
1 191 displayed, 
+3902 22 
+1455 34 
ris rig A Observe that If the number of date pairs Is W, the 
eas bed number of transform palrs 1s N/2+!. 
21989 259 
° ° 


DO YOU WISH TO EDIT ANY DATA PAIR (Y/NDE PY 


WHICH PAIR? 715 

XC 15) = 713 Date may be corrected by changing it as shown 
™ ¥C 15) = 70 to the left. 

DO YOU WISH TO EDIT ANY DATA PAIR (Y/N)! TN 


SORTING DATA... First the data Is sorted to put it In escending 
X-order. Then interpolation is performed to en- 
sure that the deta (as seen by the FFT) ts equel- 


DO YOU WISH TO SAVE THE CURRENT spaced In X. This Is a slow step. 

DATA ON DISK (Y/N)? PN If we chose to seve the data on disk (or cassette), 
we would be prompted for # file name. In the 

DO YOU WISH TO REMOVE THE cassette versions a name Is requested, but actually 

DC COMPONENT (Y/N)S ry not used, 


PLOT DATA (Y/N)? TY 


WHAT IS THE PLOT WIDTH! ?30 in most versions this refers to the number of 
hortzontal print positions. 
DISPLAY PLOT VALUES (Y/N)! 7Y You may opt to have the data pair values printed 


to the right of the plot (not available for Atari). 


DO YOU WISH AN AMPLITUDE 
VERSUS FREQUENCY PLOT (Y/N)! TY 


WHAT IS THE PLOT WIDTH? 730 


DISPLAY PLOT VALUES (Y/N)! TY 


1 
10 ° 
t 
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SSRRRSAARATAAS SERA TASS TE t .06 1 
1 t 
t 1 
eeeaseaes 2 413 a 
t : 
t 7 
aan 419 15 
t 1 
t t 
eaeeeee 3 4.25 a 
t 3 . 
: : 
as t 38 ot 
t 1 
1 + 
aan  .38 +16 
, t 
t : 
288 t 444 2 
1 t 
t 
$ 5 ° 
1 


BO YOU WISH A PHASE PLOT (Y/N)? TY 
WHAT IS THE PLOT WIDTH! 720 
DISPLAY PLOT VALUES (Y/N)? 7Y 


+ 0 +06 
t 
eeareeaeses t 406 191.25 
t : 
* t 13 22.5 
+ t 
oi t Ag 33.75 
t t 
eeeseasaesses 2 225 225 
’ + 
sax to. 546.25 
: 1 
aaen t .38 67.5 
t + 
SRSRAARAAATARS TS + 444 258.75 
, t 
* t <F 06 
I--~-I----I----I----1 


DO YOU WISH TO SAVE THESE RESULTS (Y/N)? 7H 


DO YOU WISH TG REPEAT THE ANALYSIS 
USING THE ORIGINAL DATA (Y/N)? TY 


This plot Is actually # discrete form of 

the "sinc function, but missing the 0C 

term since we chose to remove It from the 
data to start with. A ber graph pattern Is 
employed to remind you thet that the transform 
is discontinous. 


The calculation and plot ends at the Nyquist 
frequency associated with the span of the 
data (X ry ordered) and the number of deta 
points. It Is the frequency assoclate 
with the Interpolated sampling frequency. 


CURRENT DATA FILE? 


1 xa yop 
1 1 ° 
2 2 ° 
3 3 ° 
4 4 ° 
3 s ° 
é é 1 
7? 7 1 
8 8 1 
9 9 1 
10 10 1 
CONTINUE TY 

waa 1 
12 12 ° 
130013 ° 
1414 ° 
1518 ° 
16 16 ° 


DO YOU WISH TO ADD TO DATA (Y/N)! TY 
SORRY. ARRAY SPACE FILLED. 
DO YOU WISH TO DELETE DATA (Y/N)! TY 


SUGGESTION - DELETE STARTING AT THE 
END OF THE FILE. 


WHICH DATA PAIR! 716 
DO YOU WISH TO DELETE DATA (Y/N)! TN 
ERROR! THE CURRENT NUMBER OF DATA 
PAIRS IS NOT AN EVEN POWER OF 2. 
DO YOU WISH TO ADD TO DATA (Y/W)? TY 
YOU MAY ADD 1 POINTS. 
HOW MAY POINTS DO YOU WISH TO ADD! 71 
XC 16) = 716 

YC 16) = 70 
DO YOU WISH TO DELETE DATA (Y/N)! 7H 


DO YOU WISH TO EDIT ANY DATA PAIR (Y/N)? TY 
WHICH PAIR! 77. 


XC 7) = 7768 
YO?) = mh 


DO YOU WISH TO EDIT ANY DATA PAIR (Y/N)? TH 
SORTING DATA... 

DO YOU WISH TO SAVE THE CURRENT 

DATA ON DISK (Y/N) 7N 


BO YOU WISH TO REMOVE THE 
DC COMPONENT (Y/N)! TH 


PLOT DATA (Y/N)t TY 
WHAT IS THE PLOT WIDTH! ?30 
UISPLAY PLOT VALUES (Y/N)? TY. 


We are going to repeat the analysis, but make 
a change in the Input data set to test the 
Interpolation procedure. 


We can not add more date because the array has been 
filled. However, If we hed answered the original 

“number of table values" question with 32, and had 
then load a data File from disk which con 
only 16 points, more data could be added at this point. 


When @ dete pair is deleted, the pairs above shift 
down one to fill the gap. Thus the "current data 
File" table shown above becomes Incorrect above the 
patr deleted. If date pafr 7 Is deleted, pair 8 
becomes palr 7, palr 9 becomes pair 8, and so on. 
In this case we deleted the last pair in the table 
as an example, and we will later replace it. 


The program just checked the number of data pairs 
and found only 15. That Is an Invalid number for 
entry Into the FFT. You will now be given @ chance 
to correct that. 


15 points currently exist, The array size Initially 
chosen was 16. 
In this case you did not have much choice. 


We put the 16°" pair back. 


Originally pair 7 wes (7,1). We will shift that 
deta point to (7.8,1). Note thet it Is now right 
next to point B which Is at (8,1). 


The last time we removed the DC component 
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FOURIER TRANSFORM UNDER WAY.... 


FREQUENCY AMPLITUDE PHASE 
CYCLES/UNIT SCALED DEGREES 


° 1 ° 
60625 +7893 191 
+125 +308 22 
+1875 «1148 34 
2s +2357 225 
23125 +0767 36 
6375, ©1276 67 
24375 0157 239 
7 ° ° 


> wn 


ow 


10 
u 
12 
13 
4 
135 
16 


On to another example. 


10 


ee 0 80 oO 


coo eo 


Point 7 has moved so close to point 8 thet 
It Is coincident on the plot. Thus only 
point 8 Is displayed. Note thet we are stil! 
describing the same square pulse, but this 
time the points are not evenly spaced. 


Since we included the OC term this time, we can 
not meke an immediate comparison between the 
amplitude values calculated for the eventy spaced 
date and the last set. The amplitude results ere 
always scaled such thet the largest value Is 
normalized. However, observe that the phase values 
are Identical. 


il 


RUN 


In this example we will see how a dats file may be loaded from 
disk (or cassette) and more data added to It. This particular 
example Is interesting In that the data is not artificial, but 
rather the result of an experiment. tn this experiment a filter 
was probed with a train of pulses and the output displayed on 
an oscilloscope. The distorted output represents the "impulse 
response'', and its transform the transfer function. 
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THIS PROGRAM TRANSFORMS A GIVEN 
REPETITIVE DATA SET INTO ITS 
FREQUENCY DOMAIN REPRESENTATION 
USING THE FAST FOURIER TRANSFORM. 


THE USER INPUTS THE NUMBER 

OF DATA POINTS TO BE TREATED. 
THIS NUMBER MUST BE AN EVEN POWER 
OF TWO. FOR EXAMPLE? N=2°K. 

THE PROGRAM WILL THEN PLOT 

THE DATA AND ITS SPECTRUM IN 
TERNS OF AMPLITUDE AND PHASE. 


YOU MAY ALSO CHOOSE TO RECALL 
PREVIOUS DATA OR RESULTS. 


DO YOU WISH TO PLOT 
PREVIOUS RESULTS (Y/N)? TN 


INPUT NUMBER OF TABLE VALUES! 764 The actual set of data consists of 40 deta pairs. 
a The next number higher than 40 which will be accepted 
IS THE DATA TO BE LOADED Is 64, 
FROM DISK (Y/N)3 TY 
WHAT IS THE FILE NAME! TRDAR_TA Typo's can usually be corrected before the "RETURN 
are or "ENTER" key Is pressed by using "delete","rubout", 
IS THE DISK READYTY “back-arrow', or whatever the corresponding key Is 
‘on your keyboerd. 
40 POINTS WERE FOUND. Only 40 points were found. 
CAUTION! THE FILE CONTAINED The computer quickly becomes ewere thet the number 
A NUMBER OF DATA PAIRS WHICH IS of current data pairs Is not acceptable. 


NOT AN EVEN POWER OF 2. 


CURRENT DATA FILE: 


9 00018 2.8 
. 2 262 
CONTINUE ?7Y 


a1 00022 1.5 
12 100024 6 
13 100026 ° 


14 00028 716 
15 +0003 “1.2 
16 00032 “1.6 
17 00034 “1.8 
18 00036 “19 
19 00038 “ee 
20 0004 “1.7 


CONTINUETY. 


21 00042, -1.,.9 
22 00044 -1.8 
23000468 -1.6 
24 = 00048) -1..5 
25 © .0005 “1.4 
26 = 00052, -1..3 
27.0054 -1.2 
28 «00056 -141 
29° 00058 -4 
30.0006 9 

CONTINUE?Y. 

31 00062, -. 85 
32. .00064 - 8 
33.008 -t 
34.0015 3 
35.002 4 
36.0025 +35 
37.003 +25 
38.0035 2 
39.004 1S 
40,0045 at 
CONTINUETY. 


DO YOU WISH TO ADD TO DATA (Y/N)? TY 


YOU MAY ADD 24 POINTS. 
HOW MAY POINTS DO YOU WISH TO ADD? 724 


X€ 41) = 7.00003 
YC 41) 


13 


XC 42) = 7 
YC 42) = 71.06.5 


MC 43) = Ty 7 
43) = 71.95 


X( 44) = 7.00009 
YC 44) = 72.7 


This part of the print-out Is boring. 


DO YOU WISH TO DELETE DATA (Y/N)? 7H 


12 


Observe that for long time lengths and low amplitudes, 
the sampling Interval has been Increased. 


Two options exist. We can add data to get to 64 pairs, 
or we can delete data to get to the 32 level. The 
first option ts chosen. 


24 “new data points are added. They ere just the 
averages of the first 25\polnts; we are inserting 
linearly Interpolated data. 


BO YOU WISH TO EDIT ANY DATA PAIR (Y/N)I TN 


SORTING DATA... 

DO YOU WISH TO SAVE THE CURRENT 
DATA ON DISK (Y/N)? TY 

WHAT IS THE FILE NAME! ?RDATA 
13 THE DISK READY7Y 


DO YOU WISH TO REMOVE THE 
BC COMPONENT (Y/N)! TN 


PLOT DATA (Y/N)3 PY 
WHAT IS THE PLOT WIDTH! 765 
DISPLAY PLOT VALUES (Y/N): ?Y_ 


The sort takes some time for large date sets. 
remember that the sort fs followed by Interpolation 
to create an equally spaced ‘deta table, 


Recall that the file name is always required, but not 
actually used for the cassette versions. 


Always try to use @ plot width greater than the number 
of deta pairs. 
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ag graphics provided by HARMONIC ANALYZER. Since the date hes + 
as both positive and negative components, the program located : 
: the zero position and scaled the dete such that @ reasonably } 
tee baleneed graph resulted. Also, the points for short time t 
eg Intervals are very close together, end ere doubled up on : 
es print lines. When there Is doubling-up, the values are not; 
: a8 displayed to the right. ’ 
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CONTINUETY 


00064 


0015 


+002 


0025 


003 


+0035 


+0045 


35 


25 


+2 


15 


14 


FOURIER TRANSFORM UNDER WAY.... 


FREQUENCY AMPLITUDE PHASE 
CYCLES/UNIT SCALED DEGREES 


° +5604 ° 
219.7268 9953 123 
439.4535 1 “4 
459.1802. 9022" to 
878.907 8496 348 
1098.6338 .949 329 
1318.3606 .8086 312 
1538.0873  .7527 297 
1757.8141 1754 292 
1977.5408 .8377 272 
CONTINUETY, 
2197.2676 «7869 248 
2414.9944 .578 232 
2636.721 4962 225 
2856.4478 4439 218 
3076.1746 4023 209 
3295.9014 .373 200 
3515.628 .3398 188 
3735.3548 2605 168 
3955.0816 .1748 191 
4174.8082 .2069 188 
CONTINUE?TY 
4394.535 62198 163 
4614,2618 1533 152 
4833.9885 .1081 152 
5053.7153  .0896 162 
5273.4421 .0907 172 
5493.1688 41055 170 
5712.8955 41212 145 
5932.6223 .0419 113 
6152.3491  .0409 222 
6372.0757 «1 182 
CONTINUETY. 
6591.8025 0726 161 
6811.5293  .0506 161 
7031.2559 .0431 180 


DO YOU WISH AN AMPLITUDE 
VERSUS FREQUENCY PLOT (Y/N)! TY 


WHAT IS THE PLOT WIDTH! 740 


DISPLAY PLOT VALUES (Y/N)? ™ 


There are 33 transform pairs; a plot width of 40 
Is sufficient. 
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° 
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scensdessaesseusésseuses 
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saanaaeaeesaces 

daaenassaanees 

saxeuanennae 

sensdaeeene 

saxesannea 

sassenan 

anene 

auaase 


t 
waeeaee 
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tok tak - tak -Soka atta 
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I----1----1----1----I----1----1----1---- 


CONTINUETY. 


This Filter ectuelly hed no oC 
response. Host likely there was 

a small offset in the scope trace 
from which the data was obtained. 
Knowing thet there was no DC response 
we could have chosen to arbitrarily 
remove the DC component before the 
FFT calculation. That option was not 
chosen, If It had been, the amplitude 
would have been zero at zero hertz. 
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s 
t 
2 
3 
t 
t 
s 
1 
t 
t 
t 
1 
? 
t 
t 
t 
t 
3 
t 
t 
t 
t 
t 
: 
t 
1 
1 
t 
| 


° 
219.73 
439.45 
659.18 
978.91 
1098.63 
1318.36 
1538.09 
1757.81 
1977.54 
2197.27 
2416.99 
2636.72 
2856.45 
3076.17 
3295.9 
3515.62 
3735.35 
3955.08 
4174.8 
4394.53 
4614.26 
4833.98 
5053.71 
5273.44 
5493.16 
5712.89 
5932.62 
6152.34 
6372.07 
6591.98 
6811.52 
7031.25 
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DO YOU WISH A PHASE PLOT (Y/NDt TY. 
WHaT 18 THE PLOT WIDTH! 740 
DISPLAY PLOT VALUES (Y/N)? TY 


° 

I---- 

. ° ° 

: 

eegeeeneses 219.73 122.61 
t 

wees 439.45 44.29 
: 

a7 659.18 9.94 


t 
SRRAASTASAAARSTEAA AA KA ARSAAETe 878.91 347.55 


aauebantdiassesedescciaessias 1098.63 329.48 
sdavasvasssdeaseesdecesssees’ 1318.36 311.98 
saxannnenasesssoannnsnesese 1538.09 297.43 
saabaacdaveisieasaaiasaces 1737.81 292.48 
seanantsevasssiséssaseest 1977.34 271.58 
suaxanenanancessszanae 2197.27 247.99 
sansanaenneeessesanan 2416.99 231.97 
nesscenatenesstetes 2636.72 225.186 
saasientcusesstauaie 2856.45 217.67 


Saeeeeeeesasaganaee 3076.17 209.16 
t 


Ma 00 te 0 we 00 we 00 00 08 U0 te 00 Ob 00 on v0 20 wn OO 0 we oe we te te tO ene we uF be te bn On ee 20 On 00 Oo oe On ween tO C0 OF tO te OF OO On Oso teas On en ee ce te ee co on en oe 


seaeaeaereeeeeeeass 3295.9 200.19 
seassnnennsessess 3515.62 198.13 
sasseensasesess 3735.35 167.96 
qaniaceaan Baas 3955.08 181.09 
seseseasssseeeess 4174.8 188.13 
saesessnnesess 4394.53 162.86 
sacseasaeaase 4614.26 151.83 
ssssesssenenes 4833.98 152.38 
sseseeseaeaaes 5053.71 162.06 
sseseseanacesss 5273.44 172.49 
seesesnennennee 5493.16 169.59 
sensusensnese 5712.89 144.58 
sannananee 5932.62 113.36 
sanssessessesseeress 6152.34 221.93 
ertrtrrertitity 6372.07 182.43 
aesesssnnensse 6591.8 161.04 
seesnestensene 6811.52 160.76 
stesessneenenes 7031.25 180 
I----1----1---: w--1----1---- 


DO YOU WISH TO SAVE THESE RESULTS (Y/N): 7 


end of examples... 


